Quantum Discord of an Arbitrary State of Two Qubits 
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Calculation of the quantum discord requires to find the minimum of the quantum conditional 
entropy S{p ^|{n^}) over all measurements on the subsystem B. In this paper, we provide a 
simple relation for the conditional entropy as the difference of two Shannon entropies. The relation 
is suitable for calculation of the quantum discord in the sense that it can by used to obtain the 
quantum discord for some classes of two-qubit states, without the need for minimization. We also 
present an analytical procedure of optimization and obtain conditions under which the quantum 
conditional entropy of a general two-qubit state is stationary. The presented relation is also used to 
find a tight upper bound on the quantum discord. 
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I. INTRODUCTION 

Entanglement is the specific feature of quantum sys- 
tems which reveals that complete information about 
parts of a composite system does not include complete 
information of the whole system [l|, y|. However, this 
is not the only weird character of a quantum system. 
For instance, the collapse of one part of a non-entangled 
bipartite quantum system after a measurement on the 
other part is another feature unique to quantum systems. 
The quantity that captures this feature is the quantum 
discord [3, |3] . The key idea of the concept of quantum 
discord is the superposition principle and the vanishing 
of discord shown to be a criterion for quantum to classical 
transition [3|]- Furthermore, quantum discord has a sim- 
ple interpretation in thermodynamics and has been used 
in analyzing the power of a quantum Maxwell's demon 
[5| . It has also been employed in the study of pure quan- 
tum states as a resource and performance of deterministic 
quantum computation with one pure qubit [6]. The au- 
thors of Refs. 7, 8] showed that non-classical correlations 
other than entanglement can be responsible for the quan- 
tum computational efficiency of deterministic quantum 
computation with one pure qubit [6| and brought this ob- 
scured correlation measure to the spot light zone. After 
this discovery, quantum discord becomes one of the most 
frequent topics of researches in the field of quantum in- 
formation theory. Indeed, quantum discord is the differ- 
ence between two classically equivalent definitions of mu- 
tual information in the quantum mechanics language. In 
mathematical sense, discord could be obtained by elim- 
inating the whole of classical correlation from the total 
correlation measured by mutual information, by means 
of the most destructive measurement on the one party 
of the system. Mutual information of a bipartite system 
can be written as 



I{p^^) = S{p^) + S{p^)-Sip''^) 



(1) 



where p and p^ refer to the reduced density matri- 
ces of the subsystems A and B, respectively, p"^^ is 
the density matrix of the system as the whole, and 
S{p) = — Tr(p log2 p) is the Von Neumann entropy. The 
classical correlation between the parts of a bipartite sys- 
tem can be obtained by use of the measurement-base 
conditional density operator and can be written as [3| 



Cs(p"") 



sup {Sip^) 
{nf} 



5(p^|{nf})}. 



(2) 



Here the maximum is taken over all projective measure- 
ment {Hf } on the subsystem B [3], and S{p^\{U^}) = 
'YlikPkS{p^) is the conditional entropy of the subsys- 
tem A, with p^ = TrB{{l^ ® Hf ) p^^ (1-4 Hf ))/pfc 
as the post-measurement state of the subsystem A and 
Pf, = Tr(p'^^ (I'^(8)Hf )) being the probability of the A;-th 
outcome. The maximum performed in Eq. ^ can be 
taken also over all the positive operator valued measures 
(POVM) Q and these two definitions give in general in- 
equivalent results. Accordingly, discord can be calculated 
as follows 



AB\ 



„AB\ 



Db{p^^)^I{p^^)-Cb{p^^) 



(3) 
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However, one can swap the role of the subsystems A and 
B to obtain Da{p^^), which is not equal to Db{p^^) 
in general. In this paper we only consider Db{p ) and 
hence ignore the subscript B in the following. 

The optimization procedure involved in the calculation 
of quantum discord prevents one to write an analytical 
expression for quantum discord even for simple two-qubit 
systems. Quantum discord is analytically computed only 
for a few families of states including the Bell-diagonal 
states [9l .lig| . two-qubit X states [ll|,[i3|, two-qubit rank- 
2 states [ij| , a class of rank-2 states of 4 (g) 2 systems [lj| , 
and Gaussian states of the continuous variable systems 
[la |. Moreover, based on the optimization of the con- 
ditional entropy, an algorithm to calculate the quantum 
discord of the two-qubit states is presented in [iGj . It is 
also important to have some computable bounds on the 
quantum discord and some authors have obtained such 
bounds \n,\M- 



In this paper, we consider two-qubit states and ob- 
tain a simple relation for the quantum conditional en- 
tropy S{p"^^\{II^}) as the difference of two Shannon en- 
tropies. This form of the conditional entropy enables one 
to calculate its minimum for some classes of two-qubit 
states, without any minimization procedure. An ana- 
lytical progress in the minimization of the conditional 
entropy of a general two-qubit state is also presented. 
Moreover, we obtain a tight upper bound on the quan- 
tum discord of an arbitrary two-qubit state and show 
that such obtained upper bound satisfy the stationary 
conditions of the conditional entropy, when we restrict 
the measurement direction to a special subspace of M"^ . 

The paper is organized as follows. In section II, we 
consider a general two-qubit system and present a sim- 
ple relation for the quantum conditional entropy. In this 
section, we also evaluate the quantum discord for some 
classes of states, and present a tight upper bound on the 
discord of a general two-qubit state. Section III is de- 
voted to the optimization procedure. An analytical con- 
ditions under which the conditional entropy is stationary 
is presented in this section. The paper is concluded in 
section IV with a brief discussion. 



II. CONDITIONAL ENTROPY AND 
QUANTUM DISCORD OF TWO-QUBIT STATES 

A general two-qubit state can be written in the Hilbert- 
Schmidt representation as 



r,AB 



1(8)1- 



■ I (8) y • cr + y^ Ujaj ig) a^ 



(4) 

Here I stands for the identity operator, {(Ti}f^i are the 
standard Pauli matrices, x and y are coherence vectors 
of the subsystems A and S, respectively, and T = (tij) is 
the correlation matrix. Therefore, to each state p^^ we 
associate the triple {x, y, T}. Since quantum correlations 
are invariant under local unitary transformation, i.e. un- 
der transformations of the form {Ui (g) U2)p^^{Ul (g C/J) 
with C/i, U2 G SU{2), we can, without loss of generality, 
restrict our considerations to some representative class of 
the states described by less number of parameters [l9| . 
Under such transformations, the triple {x, y, T} trans- 
forms as 



Oix, y-^02y, T^OiTOl 



(5) 



where O^'s corresponds to C/i's via Ui{a-(j)U} = [Oio) -a 
[l9|. In view of this, any state of the two-qubit system 
can be written as {Ui (g) U2)p'^^{Ul ® Ul), where p^^ be- 
longs to the representative class. In the following we will 
consider a representative class such that T is diagonal, 
namely T — dia,g{ti,t2Tt^}- Concerning this represen- 
tative class, a general state of two-qubit system can be 
parameterized by nine real parameters x = (a;i,X2,a;3)*, 



V = (2/1, 2/2, 2/3)*, and T = diag{ii,i2,i3}, where t de- 
notes transposition. Accordingly, in the computational 
basis {|00),|01),|10),|11)}, a general state p^^ of this 
representative class takes the following form 



^AB 



where 



/ P\\ y\ - iy2 xx - 1x2 t\ - t2 

2/1 + iyi P22 h +t2 xi- 1x2 

Xi + iX2 h + t2 P33 2/1 - «2/2 

\ ti-t2 Xi+ iX2 2/1 + «2/2 P44 



(6) 



pil = 1 + X3 + 2/3 + ^3, P22 = 1 + 2:3 - 2/3 - ^3, 
P33 = 1 - 2^3 + 2/3 - ^3, P44 = 1 - a;3 - 2/3 + ^3- 

Now let us turn our attention on the von Neumann mea- 
surement on the subsystem B. A general such measure- 
ment can be written as 



nf = c/|fc)(fc|t/t. 



(7) 



where {|^)(A:|}^^o is the von Neumann projection op- 
erators in the standard basis of the subsystem B, and 
U e SU{2). An arbitrary element of SU{2) can be fac- 
torized as [2011 



u = n2ni, 

with Q2 and Qi defined by 



(8) 



r!2 = 



"^ 2 
e"^ sin \ 



, r!i = 



e-*''/2 



(9) 
for < 6* < vr, < < 27r, and < 7; < TT. There- 
for, we get nf = fl2\k){k\fl2 = \cr ■ nk){(T ■ nk\ where 
fiQ — —fii = h, with n — {sm9 cos(j),sm9 s[ii(f),cos9Y . 
This, explicitly, shows that only two independent param- 
eters 9 and </> are needed to characterize a general von 
Neumann measurement on the two-qubit systems. We 
can also write these orthogonal projections in the Bloch 
representation as 

nf = i(I + nfc-a), fc = 0,l. (10) 

Therefore uq and fii are coherence vectors of 11^ and Ilf , 
respectively. For further use, we calculate the expression 
nf cr^nf for fc = 0, 1 and j = 1, 2, 3, i.e. 

nf CTjIlf = |cr • nk){a ■ hk\<7j\(T ■ nk){(J ■ nk\ 

= Tr(nf a,)nf = (nfe),nf . (ii) 

Now we are in the position to perform the von Neumann 
measurement {n^}J,^Q on the qubit B. This transforms 
the total state p"^^ to the ensemble {pfc,/0^'^}J.^Q such 
that 

pr = -(i®nf)/^(i®nf), (12) 

Pk 

with Pk = Tr[(I(g)nf)p^^(I(gnf)]. By using Eqs. Q 
and (|10p in (|12l) and invoking relation pT|) we get 



p^^P^^nf, 



(13) 



where 



Pk = ^\'^ + Xk-a 



(14) 



is the post-measurement state of the qubit A, associated 
to the measurement resuh k with the corresponding prob- 
abihty 



Pk = 2 (i + y'^-fc) 



(15) 



In Eq. (|14p. the post- measurement coherence vector Xk 
is defined by 



Xk = 



X + Thk 
1 + y*-nk 



(16) 



The quantum conditional entropy with respect to the 
above measurement is defined by 

^(/|{nf})=Po^(p^|{nf})-Hpi5(pf|{nf}). (i7) 

Using 

i(l±|4|) =^(2pfe±|x-|-rnfc|), (18) 

as the eigenvalues of p^, for k — 1,2, and after some 
calculations we obtain the following relation for the con- 
ditional entropy 



s{p^\{n'^})^h,{w)-h2{po). 



(19) 



Above, and throughout this paper, h2{x) denotes the bi- 
nary Shannon entropy |21| defined by 

h2{x) = — xlog2 a; — (1 — a;)log2(l — x), (20) 

and hrn{qi, ■ ■ ■ , Qm) = - Z]"=i 1i log2 Qi is the Shannon 
entropy of the probabilities {qi, • • • ,qm}- In particular, 
hii^w) = — X)i=i ''^i lo§2 ^i is the Shannon entropy of the 
probabilities 



Wio = 



2po±\x + Tn\ 



W3A 



2pi ± \x — Tn\ 



(21) 



Note that under the transformation h — >■ —n, correspond- 
ing to 9 — >■ TT — 6 and (f) — >■ ^ ± tt, the probabilities 
([T5|) and ([2T|) transform as po ^ pi, wi i-^ w^ and 
W2 -^ it;4, leaving therefore the conditional entropy in- 
variant. On the other hand, if we perform local unitary 
transformation ^ on the density matrix, the probabil- 
ities p^ and (|2T|) . and hence the conditional entropy 
do not change provided we perform the transformation 
n — > 02n- This implies that if we find h* as the optimal 
measurement for a given state p of the representa- 
tive class, one can obtain the optimal one for any state 
p^B ^ (^u^ ^ jj^-^pAB^jj^^ ^ [/t^^ jyg^ ^y ^i^g transforma- 
tion h* -^ 02n*. 

Now the aim is to minimize the above conditional en- 
tropy. Before we give a general procedure for optimiza- 
tion, we give below some special classes of states for which 



the quantum discord can be evaluated without the need 
for any optimization. 

Bell-diagonal states. — As the first example, let us con- 
sider the Bell-diagonal states, i.e. states with the maxi- 
mal marginal entropies. For this three-parameters class 
of states we have x = 0, y = and T — diagjii, ^2,^3}- 
In this particular case we have 

which leads to the following simple form for the condi- 
tional entropy ([TI 



S{p^\{n}) = h2 



1 + \Tn\ 



(23) 



Clearly, the minimum of the above equation occurs for 
max I Tn I = max^/?^*T*^f^, which is achieved when n is 
an eigenvector of T'^T corresponding to the largest eigen- 
value t^„^ 



mmSip^\{h})^h2 



l+tr: 



(24) 



where we have defined 

imax = max{|ii|,|i2Ui3|}. (25) 

Therefore, the discord is obtained as follows 



S{p ) ^ I - hiipi, P2, P3. Pi) + h2 



where {pi}f^i are eigenvalues of p given by 

pi,2 = - (1 ± ii ± t2 - ^3) , P3A = - (1 ± ii T ^2 + <3) • 
This is in agreement with the result obtained by Luo in 

States with T^x — and y — 0. — For the second class 
we consider states such that y = and x belongs to the 
kernel of T*, i.e. T*x = 0. For this class of states we get 



\x + Th\ ^\x~ Th\ = Va;2 + h^TTh, (26) 

and therefore 
Po=Pi = 2' wi^2^W3^4 = -{l±\x + Tn\). (27) 



This leads to 



S{p^\{n}) = h2 



l + \x + Th\ 



(28) 



The minimum of the above equation occurs whenever 
|a;-fTn| = y/x^ + fi^T^Tfi takes its maximum value. This 
happens when fi is an eigenvector of T*T corresponding 
to the largest eigenvalue i^ax' therefore 



min5(p !{"•}) = /12 



l + ^x^+ tl^ 



(29) 



where tmax is defined by Eq. ((25|) . Note that the condi- 
tion T^x — requires that tiXi — for i = 1, 2, 3. Hence 
if we take, without loss of generahty, \ti\ > |t2| > l^sl > 
and set y = 0, then the states corresponding to this class 
can be obtained from the general form of Eq. ([6]) as: 
(i) States with xi = X2 = xs = 0. This corresponds 
to the Bell-diagonal states where we have studied them, 
separately. 

(ii) States with xi = X2 = t^ — 0. This corresponds to 
a three-parameter subclass of the so-called X states. In 
this case, the discord is obtained as 



S{p ) ^l~hi{fJ.l,^l2,^J-3,^J-4:) + h2 



1 + V4 + '' 



and 



Po=Pi = ^, wi^2 ^W3^4^ -{l±\x + Teo\) . (35) 



Therefore 



min ^(/|{n}) < min^(p^|{eo}) 



= min /i2 



l + y^x^+ e*Tr*eo 



l + v/^^+i2 



(36) 



where {/ii}^^x ^'"^ eigenvalues of p given by 

1 

4 
1 



Mi,2 = 7 1 i±x/{ti + t2y + xi,, 



mA^^\'^±\/iti-t2y 



(iii) States with xi — t2 — ts — 0. This corresponds to a 
three-parameter subclass of the zero-discord states, 
(iv) States with ii = ^2 = ^3 = 0. This also corresponds 
to a three-parameter subclass of the zero-discord states. 



where t^ — maxegT'Teo with the maximum taken over 
all unit vectors eo G TZ^. Evidently, if two vectors 
T*x and y be nonzero and linearly independent then 
dim??.-'- = 1, so that the unit vector gq € TZ^ is unique. 
Otherwise dim 7?,^ > 1, so we can choose cq G TZ^ such 
that Iq — max e^T'^Teo, giving a tighter bound. This 

completes the proof of ^. Using Eq. ^ in Eqs. @ 
and ([3]), we obtain the desired bounds ([5^ and (1551) . ^^~ 
spectively. 



Interestingly, the above two examples motivate us to 
introduce an upper bound for the quantum discord. Let 
TZ be the subspace of R^ spanned by T^x and y, i.e. 
TZ = span{r*x, y}, and let TZ^ denotes the orthogonal 
complement of TZ, i.e. the set of all vectors in R"^ that 
are orthogonal to every element of TZ. Hence we have 
TZ + TZ^ =R'\ 



Theorem 1 The conditional entropy I119\) is bounded 
from above as 



min^(/|{nf})</.2(ii#±^K (30) 



{nf} 
where x — \x\, and 



tl = max eiT^Teo- 



(31) 



Accordingly, the classical correlation and the quantum 
discord have the following lower and upper bounds, re- 
spectively 



C{p^")>S{p^)-h2 



l + ^x^Ttl 



(32) 



Q{P^^) < Sip^) - S{p^^) + h2 ( i±4^±l 1(33) 



Remarkably, the above bound is tight in the sense that 
for all states that TZ-^ = R'^, the bound is achieved. This 
happens, for example, when x = y = OoT:T*x = y = 0. 
In these cases tg becomes the largest eigenvalue of r*T 
and eo is the corresponding eigenvector. As we will show 
in the example below, it may happens TZ-^ ^ M? but 
the upper bound pO]) is achieved. In these cases the 
absolute minimum of the conditional entropy happens 
for directions in TZ^ C R-^. Latter, we show that the 
above upper bound corresponds to the stationary point 
of the conditional entropy when we restrict ourselves to 
the solutions in TZ^ . 

Recently an upper bound for the quantum discord is 
obtained in [13 as Q{p^^) < S{p^). A comparison of 
this with the upper bound presented in (j33p shows that 

for all states for which S{p^^) — /i2 



i±4^) >o. 



our bound is stronger. 

As an illustrative example let us consider a two- 
parameter class of states discussed in [22] 



^^>'^)-2l 
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-a + b 







a 



, (37) 



Proof Let us concern about all measurement directions 
eo living in TZ^, i.e. j/*eo = 0, (T*f)*eo = 0; then 



where < a < 1 and a— l<b<l — a. The discord of 
this state is [221 



\x + Teo\ = \x- Teo\ = ^x^ + e^TTeo, (34) 



g(p'^^(a,6)) =min{a,Q}, 



(38) 



where 
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For this state we get 



-y 




T ■ 



a 
-a 
2a - 1 



(39) 



(40) 



Clearly, T*x = —(2a — l)y, so that for all values of a 
and 6 y^ we have TZ — spanly}. Therefore an arbitrary 
element of TZ^ can be written as eo = (cos 0, sine/), 0)*. 
In this case we get egT*Teo = a?, which is independent 
of (j), so that to = o? ■ This means that every vector in 
the two-dimensional subspace Ti-^ , corresponding to the 
xy-plane, gives the desired upper bound. Therefore, we 
obtain 



and 



min S{p |{eo}) = /i2 



S{p^) = 5(p^) 



'l + Va2 + 62^ 



5(p^«) = /^3 «, 



1 — a — 6 1 — a + 6 



(41) 

(42) 
(43) 



where /i3(mi, ^^2, A*3) = - Y.i=i f^-i log2 Mi- After some cal- 
culations wc find that the inequality (|33p leads to 



g(p^^(a,6))<«, 



(44) 



where q is defined in Eq. ([M]). A comparison of this 
with Eq. ((38|) shows that for some region of parameters, 
namely when q < a, our upper bound is tight. Figs. [T] 
and [2] illustrate this fact. These figures also reveal that 
for this class of states our bound is stronger than the 
upper bound of Ref. [13 ■ 



III. OPTIMIZATION 

In this section we present an analytical procedure for 
optimization of the conditional entropy for a general two- 
qubit state. We also provide examples for which one 
can obtain the the minimum, analytically. In order to 
determine the minimum of the conditional entropy ([T9| , 
we have to calculate its derivatives with respect to 9 and 
(/>. To do this we need to calculate derivatives of the 
probabilities given by Eqs. ([T5|) and (PT|) with respect to 





















as 

1 
1 "■'' 






1 

J »' 

|04 

Q 








^-r" 




3.0 0,2 0.4 0. 


O.S 






1.0 ai 0.2 0.1 


04 


5 0.6 






















1 "'' 

1.4 
Q 
□ .2 


M 




O.S 

1 '■" 

|o4 






\ 


^^^^^^^^^ 







.0 0.1 O.Z OJ 0.4 0.5 


00 0.02 0.04 


0.06 


0.0! 0.1 



FIG. 1. (Color online) Quantum discord [dashed-blue lines], 
our upper bound [solid-black lines], and the upper bound of 
Ref. [i3] [dotted-red lines] are plotted versus a for the state 
p^^{a,b) with: (a) h = 0.1, (b) 6 = 0.3, (c) 6 = 0.5, and (d) 
6 = 0.9. 
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FIG. 2. (Color online) Quantum discord [dashed-blue lines], 
our upper bound [solid-black lines], and the upper bound of 
Ref. [13] [dotted-red lines] are plotted versus h for the state 
p'^^(a,6) with: (a) a = 0.1, (b) a = 0.3, (c) a = 0.5 (d) 
a = 0.9. In the cases (c) and (d), our upper bound completely 
coincide with the quantum discord. 



and (j). For instance, their derivative with respect to 6 
are as follows 



dpQ 



1 ^l^t ^ 



with 






Th + x 

\Th + x\ ' 



y±T*Z+ 
-y±T*Z 



Th — X 
\Tfi — x\ 



(45) 
(46) 
(47) 

(48) 



and the unit vector n^g is defined by 
dn 



nft 



de 



(cos COS 0, COS sin 0, — sin 6*) . (49) 



Evidently h- h g = Q. By defining the nonunit vector fi , 

by ' ' 



dh 



",0 



(— sin6'sin0, sin0cos0, 0)*, (50) 



orthogonal to both n and hg, we get a similar equa- 
tions for the derivatives of the probabilities with re- 
spect to (j), but now fi^g is replaced by fi^^. Using 
the above equations, we arrive at the following condi- 
tions for the stationary conditions dS{p'^\{T\-^\)/d9 — 

a5(p-4|{nf})/a0 = o 



n*fl A = 0, 



n\ 



A = 0, 



(51) 



where A is a vector defined by 



log; 



W1W2PI 



y+ 



log; 



^^1 
W2 



T*Z+ + 



log; 



W3 



Conditions (ISTt can be written also as 



h\ 



A = 0, 



(53) 



where h± is any vector perpendicular to n, i.e. n^n = 
0. Equation (15^ shows that the stationary points of 
S{p^\{Il^}) happens whenever the vector A is in the 
direction of h. If we proceed the optimization process 
by using the Lagrange multiplier A, having n*h ~ 1 
as a constraint, we find for the stationary condition 
d [S{p^\{n^}) - A(n*n - 1)] = the following relation 



A = An, 
where A is defined in Eq. (j52[) . and A is given by 



(54) 



A^ 



Wi 



l0g2 



W2 



f*Z^ 



, W3 

W4 



x'Z- 



45(p^|{nf}) 



log; 



1011021113104 
PoPi 



(55) 



The stationary condition (j54p is equivalent to that given 
by Eq. (|53p in the sense that both conditions require 
that in the extremum points, the vector A should be di- 
rected to n. Unfortunately, these conditions do not have 
simple solutions, for A as well as A depends also on n. 
Moreover, knowing the extremum points of the condi- 
tional entropy is not enough to establish its minimum, 
and we need, in addition, to evaluate the conditional en- 
tropy in the extremum points. Below we exemplify these 
conditions for some particular cases where we have al- 
ready obtained the minimum of the conditional entropy 
without optimization. 

Bell-diagonal states. — For this class of states we al- 
ready showed that the optimal measurement lies in the 



direction of the eigenvector of T'^T, corresponding to the 
largest eigenvalue. In this case the vector A takes the 
following form 



A 



Th 



log2 



1 + \Tn\ 
1 - \Th\ 



T^Th. 



(56) 



It is clear that the condition ([5^ leads to n^^T^Tn — 0, 
which has solutions when h is an eigenvector of T'T. But 
the 3x3 matrix T*r has three nonnegative eigenvalues 
{if,t|,i|} corresponding to the eigenvectors {61,62,63}. 
Therefore, the conditional entropy in the stationary 
points is obtained as 



S{p%ek})^h2 



l + \tk\ 



for A: =1,2, 3. (57) 



Simple evaluation shows that the minimum of the above 
equation happens when \tk\ corresponds to the largest 



In this case the vector 





eigenvalue of VT^T. 






States with Tx = and y — 0. — In tl 




A takes the following form 


T^Z-. 

(52) 


A- 2 

\Tfi + x\ 


l + \Th + x\' 
1 — \Tn + x\ 



T^Th. 



(58) 



Again the condition ([55)) leads to fi^j^T^Th = 0, which has 
solutions when h is an eigenvector of r*r. Therefore, in 
the stationary directions {61,62,63} we get 



5(/|{e,.}) = h2 



1 + V^ 



'+tl 



for A: =1,2, 3. 



(59) 
Again, evaluation shows that the minimum of the above 
entropy happens when t^ corresponds to the largest 
eigenvalue of r*T. 

Examples above show that the stationary condition 
given in Eq. (j53|) gives us, in general, more than one 
solution for the measurement direction n, and we have 
to find the optimal one by further evaluations. 

Let us mention here that the upper bound expressed 
in the theorem [T] can be actually obtained from a special 
class of solutions of the stationary condition ([5^ . namely 
the one restricted to the subspace TZ-^ of M? . To show 
this, we have to show that if the measurement direction 
60, living in the subspace TZ^ of E'^, satisfy the stationary 
condition ([55]) then the upper bound (1501) is obtained 
as the minimum of the conditional entropy. Concerning 
about all measurement directions eo living in TZ^\ then 



5(/|{eo}) = h2 



l + \x + TeQ\ 



and 



A 



\Teo 



logi 



1 + irgp + x\ 

1 - iTeo -I- x\ 



T'Tio, 



(60) 



(61) 



where by assumption, cq should satisfy the stationary 
condition (1551) . Now if diuiTZ = 2, i.e. when two vec- 
tors T*x and y are nonzero and independent, then the 



stationary condition (|53l) can be expressed as 



IV. CONCLUSION 



if A ^ 0, [T^xfA ^ 0. 



Using (pT|) in the above conditions, we get 



y*r*Teo = 0, 



(r*x)*r*reo = o, 



(62) 



(63) 



where imphes that T*Teo € TZ-^. But since dim 7?."^ = 1, 
therefore T^Tcq should be in the direction of eo, or equiv- 
alently eo should be an eigenvector of T^T corresponding 
to the eigenvalue tg. On the other hand, if one or both of 
vectors T*x and ybe zero or they are nonzero but linearly 
dependent, then 1 < dim 7?.^ < 3. In this case the two 
conditions given by Eq. (I63p are no longer independent 
and eo is not unique. In any case we get 



min 5(p^|{eo}) 



i + V^^Ttf 



(64) 



where ig is defined by Eq. (pit . Equation ([M| gives 
the conditional entropy in the stationary point eo, re- 
stricted to the subspace TZ-^ G R"^, and is equal to the 
upper bound for the conditional entropy given by Eq. 
(l30t . Evidently, if diuiTZ = 0, or if the conditional en- 
tropy evaluated in 7?. € R'^ does not be smaller than Eq. 
()64|) . then the minimum given by Eq. (|64|) gives the ab- 
solute minimum of the conditional entropy and the upper 
bound (l30l) is achieved. 



All difhculties in calculating the quantum discord arise 
from the difficulty in finding the minimum of the quan- 
tum conditional entropy S{p^^\{Il^}) over all measure- 
ments on the subsystem B. In this work, we have pre- 
sented a simple relation for the quantum conditional en- 
tropy of a two-qubit system, as the difference of two 
Shannon entropies. Using it, we have obtained the quan- 
tum discord for a class of states for which the conditions 
r*x = and y = are satisfied. This class of states 
includes the Bell-diagonal states, a three-parameter sub- 
class of X states, and some zero-discord states. For these 
states, it is shown that the quantum conditional entropy 
reduces to the binary Shannon entropy, so that in mini- 
mizing the conditional entropy we encountered with the 
simple problem of minimizing binary Shannon entropy. 
We have shown that such obtained relation for the con- 
ditional entropy can be used to present a tight upper 
bound on the quantum discord. We have exemplified 
this bound for a two-parameter class of states and show 
that in some regions of the parameters our bound is sat- 
urated. A general procedure of optimization of the con- 
ditional entropy is also given, and by some examples it 
is shown that it can be used to calculate the quantum 
discord for some classes of states, analytically. It is also 
shown that such obtained upper bound on the conditional 
entropy is actually a solution of the stationary conditions 
of the conditional entropy, when we restrict the measure- 
ment direction to a subspace of M.^. We hope that our 
results provide further insight into the problem of quan- 
tum discord. The method presented in this paper can be 
generalized to higher dimensional systems. 
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